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ABSTRACT 

Unitarity of evolution in gravitational collapses implies existence 
of macroscopic stable horizonlcss objects. With such objects in 
mind, we study the effects of anisotropy of pressures on the sta¬ 
bility of stars. We consider stars in four or higher dimensions and 
also stars in M theory made up of (intersecting) branes. Taking 
the stars to be static, spherically symmetric and the equations of 
state to be linear, we study ‘singular solutions’ and the asymp¬ 
totic perturbations around them. Oscillatory perturbations are 
likely to imply instability. We find that non oscillatory pertur¬ 
bations, which may imply stability, are possible if an appropriate 
amount of anisotropy is present. This result suggests that it may 
be possible to have stable horizonlcss objects in four or any higher 
dimensions, and that anisotropic pressures may play a crucial role 
in ensuring their stability. 
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1. Introduction 


A sufficiently massive star in four dimensional spacetime is believed to 
be unstable against gravitational collapse. Depending on its mass and its 
evolutionary stage, the star may collapse to form a white dwarf or a neutron 
star or, if supermassive, form a black hole. White dwarfs or neutron stars 
may collapse further if they gain sufficient mass, for example, by accretion. 
Thus, all sufficiently massive objects are expected to collapse to ultimately 
form black holes [1]. The same is expected in higher dimensional spacetime 
also. 

Assume that a black hole is formed in a collapse. It has a horizon, emits 
ffawking radiation, and is believed to evolve unitarily. Presence of horizon 
means that no information from inside the horizon is accessible to an ontsicle 
observer. But unitary evolution means that inside information must become 
accessible to an outside observer atleast after some time, fforizon must then 
cease to exist from this moment of accessibility. The black hole at this 
time should still be of macroscopic size so that information density within is 
atrnost of order Planckian scale, and not parametrically larger [2, 3]. And, 
the time when the horizon ceases to exist is expected to be the Page time 
which is of the order of evaporation time when about half the black hole has 
evaporated [4, 5, 6, 7]. 

Another scenario is also possible where black holes do not form at all. 
The gravitational collapse, that would have led to a black hole, will instead 
lead to an object which has no horizon as, for example, in Mathur’s fnzz ball 
proposal [8, 9, 10]. Unitary evolution in this scenario will proceed as for any 
system with large number of degrees of freedom. 

Thus, whichever scenario proves to be correct, the unitarity of evolution 
implies that there must exist horizonless objects which form about half way 
through a black hole evaporation if a black hole has initially formed in a 
collapse, or right after a gravitational collapse in which a black hole might 
have been expected to form. 

It must then be possible to construct such horizonless objects using ap¬ 
propriate sources just as, for example, one constructs neutron stars using 
Oppcnheimer - Volkoff equations and appropriate equations of state. In the 
case of neutron stars, one has a very good knowledge about the nature and 
the properties of the constituents, namely protons, neutrons, quarks, glu¬ 
ons, et cetera. In the case of horizonless objects, which can be macroscopic, 
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such knowledge is absent at present. These horizonless objects may be the 
massive remnants suggested in [2] whose size and mass depend on the in¬ 
formation contained; or they may be the ones which develop ‘firewalls’ as 
suggested in [11, 12, 13]; or they may be the ‘fuzz balls’ of Mathur’s proposal 
which have no horizon at anytime [8, 9, 10]. Recently, we have also argued 
for the existence of such horizonless fuzz ball like objects in any quantum 
theory of gravity where singularities are resolved and evolutions are unitary 
[14, 15, 16]. A variety of horizonless objects have been advocated in the past 
also from several points of view, see [17] - [28] for a sample of them. 

Although nothing is known rigorously about the horizonless objects and 
their constituents, it is physically reasonable to expect them to have the fol¬ 
lowing properties, (i) Their central densities are likely to be of the order of 
Planckian densities since the quantum gravity effects are expected to play 
an important role in creating them, (ii) Their sizes are likely to be of the 
order of Schwarzschild radii since they are expected to be similar to black 
holes, only without horizons, (iii) Being the end products of prior gravi¬ 
tational collapses and evolutions, they must be stable against any further 
gravitational collapse, (iv) Their constituents must have a large number of 
entropic degrees of freedom and, hence, thermodynamics must be applica¬ 
ble. Therefore these constituents, even if highly quantum in nature, may be 
modelled using density, pressures, and suitable equations of state, (v) These 
horizonless objects can be arbitrarily massive. Hence, atleast qualitatively, 
the corresponding spacetime may be described by general relativity equa¬ 
tions with appropriate energy momentum tensors and constituent equations 
of state. 

In this paper, we consider static, spherically symmetric cases and focus 
on property (iii), namely the stability property. In the following, for the sake 
of brevity, we will refer to these horizonless objects as stars although they 
are not stars in a conventional sense. For example, their central densities are 
likely to be of the order of Planckian densities; and their sizes are likely to 
be of the order of Schwarzschild radii. 

Sufficiently massive conventional stars have been found to be unstable 
against gravitational collapse. But, in his study of relativistic stars in arbi¬ 
trary dimensions [29], see also [30], Chavanis found among other things that 
stars can be stable against collapse, albeit only in higher dimensions. To the 
best of our knowledge, this is the only work which has shown the possibility of 
arbitrarily massive stars being also stable; moreover, the radii of these stable 
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stars scale as Schwarzschild radii. 1 Chavanis considered static spherically 
symmetric stars made up of a perfect fluid with a linear equation of state 
and enclosed them in a box. 2 Following the methods of Chandrasekhar [34], 
he then obtained ‘singular solutions’ and analysed the asymptotic pertur¬ 
bations around them. He has shown in detail that the damped oscillatory, 
or monotnic non oscillatory, behaviour of these perturbations around singu¬ 
lar solutions leads to the damped oscillatory, or monotnic non oscillatory, 
behaviour for the mass - central density profile of the star in the asymp¬ 
totic region of large central density or large radius. He has shown further 
that equilibrium configurations become unstable beyond the first maximum 
in this profile, hence its oscillatory behaviour in the asymptotic regions will 
imply instability against collapse. Studying the effects of higher dimensions, 
he then found that the mass - central density profiles of the stars in eleven 
or higher dimensions are non oscillatory in both the asymptotic and non 
asymptotic regions, and imply stability. 

If the number eleven above counts both compact and non compact direc¬ 
tions then stability in, for example, four dimensional non compact spacetime 
may be obtained by having a eleven dimensional spacetime with seven com¬ 
pact dimensions, which is quite natural in the context of M theory. With this 
motivation, in [37], we generalised the work of Chavanis to include compact 
directions and multi component perfect fluids, the later also appearing quite 
naturally in string and M theory black branes. However, we found that, even 
in these generalised cases, eleven non compact dimensions are needed for 
stability. In particular, this implied that stars in four dimensions and those 
in M theory are unstable and that they will collapse if sufficiently massive. 

Reviewing the underlying assumptions in [29, 37] and studying carefully 
some of the works on the horizonless objects, particularly [19] - [24], we 
realised that the pressures inside the stars along the noncompact spatial 
directions need not be isotropic. The assumption about isotropy may be 
unwarranted and restrictive : it is not required by spherical symmetry and, 

1 Similar results were also obtained around the same time in [31, 32], but in asymptot¬ 
ically anti de Sitter spacetimes whose radii limit the masses of the stable stars. 

2 This box is needed since the radius of the star will be infinite otherwise [33]. As ex¬ 
plained in [29], enclosing the star within a box prevents the evaporation of its constituents 
and makes its total mass finite. The radius and the mass may also be made finite by 
constructing a composite configuration consisting of a perfect fluid core and a crust of 
constant density or a gaseous envelope exerting a constant pressure, see [34, 35, 36]. 
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in general, the pressures may be different along the radial and the transverse 
spherical directions of the non compact space. For conventional boson stars, 
anisotropic pressures arise naturally when scalar fields are present [38, 39]. 
See also the reviews [40, 41, 42, 43] and the recent paper [44], For the 
horizonless objects, which are implied by unitarity in lieu of black holes 
but are referred to as stars here, there may be other mechanisms giving 
rise to anisotropic pressures. However, since not much is known about the 
constituents of such objects, we do not know any possible exact mechanisms. 

In this paper, therefore, we simply assume that pressures in the stars 
may be anisotropic and study, a la Chavanis, the effects of such anisotropy 
on the stability of stars. Following closely the works in [29, 34, 37] and, gen¬ 
eralising them now by including anisotropic cases, we study stars in higher 
dimensional spacetime which may have compact toroidal directions also. The 
stars are assumed to be static and spherically symmetric in the non compact 
space and to have suitable isometries along the compact directions. Taking a 
suitable ansatz for the metric, we write down the equations of motion. Then, 
as in [29, 34], we assume linear equations of state and obtain the ‘singular 
solutions’ and the asymptotic perturbations around them. The oscillatory 
or non oscillatory behaviour of these perturbations lead to corresponding 
behaviour for the mass - central density profile in the asymptotic regions. 
We assume that, as in Chavanis’ works [29], the equilibrium configurations 
become unstable beyond the first maximum in the mass - central density 
profile; hence that the clamped oscillatory behaviour of this profile implies 
instability whereas monotnic non oscillatory behaviour throughout implies 
stability. Studying then the asymptotic perturbations around the ‘singular 
solutions’, we obtain the criteria under which the perturbations are non oscil¬ 
latory. These will then be the necessary criteria for stability. 3 We perform 
the above analysis first for stars in spacetime with no compact directions; 
then include compact directions; then repeat the analysis for two examples 
of stars in M theory made up of stacks of (intersecting) two branes and five 
branes. The formulation presented in this paper may, however, be used to 

3 One must now show that the equilibrium configurations in the mass - central density 
profile become unstable beyond the first maximum. And, for the case where asymptotic 
behaviour is non oscillatory, one must show that this profile remains monotonic every¬ 
where, including the non asymptotic region also. Showing these is beyond the scope of the 
present paper since it requires a knowldege of the nature and the properties of constituents 
of the horizonless objects, which is lacking at present. 
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study a variety of other examples also. 

We find that, in an (m + 2) dimensional non compact spacetime with 
m > 2 , non oscillatory perturbations are possible for any value of m if an 
appropriate amount of anisotropy is present, namely if 

fl — p 9 — m 
fl — 4m 

where II and p are the pressures along the radial and the transverse spherical 
directions. Also, we find that presence of compact directions does not change 
this result. Similar result follows for stars in M theory also. It can now be 
seen that if m > 9, namely if the non compact spacetime is eleven or higher 
dimensional, then the isotropic case (above ratio = 0) is included in the above 
range, thus reproducing the result in [29]. For lower values of m, certain 
amount of anisotropy is needed to obtain non oscillatory perturbations; for 
example, for four dimensional spacetime, m = 2 and the above ratio needs 
to be > | . 

This result suggests that it may be possible to have stable horizonless 
objects in four or any higher dimensions, and that anisotropic pressures may 
play a crucial role in ensuring their stability. Although much remains to be 
done, for example resolve the issues mentioned in footnote 3, it is worth em¬ 
phasising that this is an important result because it bears on the horizonless 
objects which are implied by unitarity in lieu of black holes, and it points out 
a necessary ingredient for their stability. To actually construct such objects, 
however, requires detailed understanding of many issues such as the nature of 
the constituents and the physical mechanisms that may provide the required 
amount of anisotropy. We will discuss briefly these and other issues at the 
end of the paper. 

This paper is organised as follows. In section 2 , we first mention briefly 
the relevant aspects of [29] which are used here. We then briefly describe M 
theory stars and mention how the equations of state are obtained. In section 
3, we set up our notations and conventions and present the equations of 
motion in a suitable form. In section 4 , we give the details of the asymptotic 
analysis, which involve singular solutions and perturbations around them, 
and outline briefly the significance of perturbations. In section 5, we obtain 
singular solutions for stars in spacetime with no compact directions, study the 
perturbations around them, and obtain the criteria for their non oscillatory 
behaviour. In section 6, we repeat the analysis for two examples of stars in 
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M theory : in one, stars are made up of a stack of two branes or five branes; 
in another, stars are made up of four stacks of intersecting branes, two stacks 
each of two branes and five branes. In section 7, we conclude with a brief 
summary and a discussion of some of the issues which require further study. 

2. Some general remarks 

In this section, we first mention briefly some aspects of [29] which are 
relevant here. For more details, see [29] and the references there. We then 
briefly describe M theory stars and mention how the equations of state are 
obtained. 

Chavanis considers static spherically symmetric star with its constituents 
obeying linear equations of state. In cosmological contexts, it is a standard 
practice to use linear equations state. In the contexts of stars, Chavanis ex¬ 
plains several situations where such equations of state arise and the reasons 
for using them. The star is then enclosed in a box of radius r* which prevents 
the evaporation of its constituents and makes its total mass finite. The gen¬ 
eral relativistic equations are then solved to obtain hydrostatic equilibrium 
configurations. 

Considering a series of equilibria, the mass - central density profile is 
obtained and shown to have damped oscillatory behaviour. The hydrostatic 
equilibrium configurations become unstable beyond the first maximum in 
this profile. This is also shown to correspond to conditions for nonlinear 
dynamic stability. Thus, for a given volume, the star has a maximum mass 
above which it becomes unstable. 

Following Chandrasekhar [34], singular solutions and the asymptotic per¬ 
turbations around them are obtained. The singular solutions give the scal¬ 
ing relations between mass, radius, and other quantities of the stars. The 
damped oscillatory, or monotonic non oscillatory, behaviour of the pertur¬ 
bations are shown to be responsible for the corresponding behaviour for the 
mass - central density profile in the asymptotic region. 

Chavanis obtains the mass - central density profile in the non asymptotic 
regions also. He shows that the profile in D < D crit dimensional space- 
time starts from origin, increases monotonically, reaches a maximum, then 
decreases, and oscillates with damped amplitudes around the constant value 
line given by the singular solutions; whereas, in D > D cr i t dimensional space- 


7 



time, the entire profile is monotonically increasing, and asymptotes to the 
constant value line given by the singular solutions. The critical dimension 
D crit ~ 11 and depends on the pressure to density ratio. See Figure 23 in 

[29]- 

If the mass - central density profile is oscillatory then the scaling be¬ 
haviour given by the singular solutions is not stable since the stability is lost 
beyond the first maximum of the profile. The singular solutions are then 
of no relevance. The mass - radius relations and other quantities of stable 
stars must be obtained from detailed analysis of the equations in the non 
asymptotic regions before the first maximum. 

If the entire mass - central density profile is monotonic and non oscillatory 
then the series of equilibria represented by its points are all stable. Then the 
singular solutions correspond to stable configurations and they can be used 
to obtain mass - radius relations and other quantities of stable stars in the 
limit of large central density or large radius. 

In the case of horizonless objects which are implied by unitarity and which 
we study here, nothing is known rigorously about the nature and properties 
of their constituents. Hence a thorough analysis as in [29] is presently not 
possible. Nevertheless, it is possible to make some progress using the phys¬ 
ically reasonable properties of such objects listed in the Introduction. We 
assume a linear equation of state for the constituents. The stiffest equation 
of state (pressure = density , sound speed = light speed ) is of this type. 
Also, this linearity helps in obtaining analytically the singular solutions and 
the asymptotic perturbations around them, and in studying whether mono¬ 
tonic non oscillatory behaviour is possible. We enclose the star in a box of 
radius r* which prevents the evaporation of its constituents and makes its 
total mass finite. Then, as in [29], the damped oscillatory, or monotonic 
non oscillatory, behaviour of these perturbations can be shown to lead to the 
corresponding behaviour for the mass - central density profile of the star in 
the asymptotic region. However, we can not obtain this profile in the non 
asymptotic regions, nor study the stability properties around the maxima, 
without knowing the detailed properties of the constituents. 

Detailed knowledge of the constituents of the horizonless objects is also 
needed in order to understand the fate of, for example, a collapsing massive 
neutron star which would have formed a black hole in the standard scenario. 
This knowledge is needed to understand what happens to neutrons, protons, 
et cetera in the quantum gravity regime where the singularities are assumed 



to be resolved. This requires understanding the relation between quantum 
gravity theory and standard model particles. 

M theory stars 

In String/M theory, the entropy and Hawking radiation of a class of ex¬ 
tremal and near extremal black holes have been understood rigorously in 
terms of various intersecting brane configurations and the low energy exci¬ 
tations on them. Mathur’s fuzz ball proposal also arises naturally in this 
context. The entropy and Hawking radiation of neutral or far from extremal 
black holes are not understood equally rigorously, but their explanations are 
likely to be in terms of intersecting brane antibrane configurations along 
the lines given in [45] - [54], In the intersecting brane configurations, var¬ 
ious stacks of branes wrap around compact toroidal directions, intersecting 
according to BPS rules whereby, in M theory, 4 two stacks of five branes 
intersect along three common spatial directions; a stack each of two branes 
and five branes intersect along one common spatial direction; and two stacks 
of two branes intersect along zero common spatial direction. See [55] and 
the references therein. Given the understanding of near extremal black hole 
properties in terms of intersecting branes, and given the fuzz ball proposal, 
it is natural to expect that there must be stars in M theory made up of in¬ 
tersecting branes 5 which, when sufficiently massive, will collapse and form 

4 The brane configurations in string and M theories are equivalent and are related by 
chains of U duality operations involving dimensional reduction and upliftment between 
string and M theory, and the T and S dualities of the string theories. Hence, in the 
following, we will restrict ourselves to M theory. The corresponding string theory results 
are straightforward to obtain. 

5 M theory stars which are considered here and are made up of, for example, two branes 
may be visualised as follows. The spacetime is eleven dimensional with two compact 
toroidal spatial directions. The non compact space is eight dimensional, described by 
polar coordinates (r, 6 \, •••, # 7 ) . Thus there is a torus at every point of the non 
compact space. Stacks of two branes wrap around these tori such that there are a total of 
n(r) number of two branes per unit volume of the non compact space. This number density 
n(r) ^ 0 for r < r* and = 0 for r > r* . This provides a picture of an M theory star of 
radius r* made up of two branes. This picture is analogous to that of a four dimensional 
star of radius r* where (r, 0 ±, 62 ) are the polar coordinates and n(r) is the particle density. 
The stack of branes wrapping the torus at a point is like a ‘particle’ in the non compact 
space. One may similarly visualise M theory stars made up of other intersecting brane 
configurations. A ‘particle’ in the non compact space now is the A f stacks of two and five 
branes wrapping the tori at a point, with necessary isometries, and intersecting according 
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either black holes having horizon or fnzz ball like objects having no horizon. 
Also, for the same entropic reasons which are explained in [56, 57, 58] where 
early universe was studied in string/M theory, we may assume that stars in 
M theory are made up of stacks of intersecting branes. 

M theory stars may be studied using the formalism given in this paper. 
Assuming the spatial directions of the brane worldvolumes to be toroidal and 
assuming necessary isometries, the M theory brane configurations consisting 
of J\f stacks of two and five branes which intersect according to BPS rules, can 
be described by A f seperately conserved energy momentum tensors Tmn(i), 
I = 1, 2, - - •, A f , with appropriate equations of state among their compo¬ 
nents. See [59, 60, 61] and the references therein. The equations of state 
may follow from an action; or they may be derived using the microscopic 
dynamics of the constituents which are far from being rigorously known; or 
they may simply be postulated as an ansatz. 6 

In M theory, equations describing black hole spacetimes follow from an 
eleven dimensional low energy effective action [45, 46, 55, 59, 60, 61]. In 
the context of an expanding universe, the equations of state for intersect¬ 
ing branes have been derived in certain approximations using microscopic 
dynamics of branes [56, 57]. In the context of an expanding universe, and 
also of stars, we used U duality symmetries and derived a relation among the 
components of the energy momentum tensor [58]. ' These relations, one each 
for each of the J\f stacks in the intersecting brane configurations, follow as 
a consequence of U duality symmetries and, therefore, must always be valid 
independent of the details of the equations of state. The equations of state 
in [45, 46, 55, 59, 60, 61, 56, 57], which were obtained by other methods, all 
obey these U duality relations. 

to BPS rules. This is a rough picture; the details of how the ‘particles’ interact are far 
from being rigorously known. 

6 This situation is similar to that for, for example, charged black holes versus for stars 
made up of charged particles or for an expanding universe containing them: In the case 
of black holes, equations are obtained from an action. In the case of stars, or an expand¬ 
ing universe, one may use statistical mechanics to describe charges, anticharges, and the 
photons between them; or one may simply take, for example, a linear equation of state as 
an ansatz. 

7 Using dimensional reduction and upliftment, T dualities, and S dualities, we obtained 
the U duality relations first in the cosmological context in [58] and then, in unpublished 
notes, for stars and black branes also. See the comments and the application of these U 
duality relations in [62]. 


10 



The U duality relations need to be supplemented with the equations of 
state in, say d dimensional, non compact part of the total eleven dimensional 
spacetime. These data on the equations of state are same as those needed 
for any expanding universes or for any stars in d dimensional non compact 
spacetimes but now with no compact dimensions. See [58, 63, 64, 65, 37] 
for examples and more details, as well as for applications and non trivial 
consequences of these U duality relations for expanding universes and for 
stars in M theory. 

Note that, in the context of M theory also, understanding the fate of, for 
example, a collapsing massive neutron star which would have formed a black 
hole in the standard scenario, requires understanding the relation between 
M theory and standard model particles. 

3. Equations of motion 

In this paper, we consider static cases which are spherically symmetric 
in higher dimensional spacetime. We will also consider such stars in eleven 
dimensional M theory. Although the formalism and the results in this paper 
are applicable to conventional stars also, our main interest here is in the 
horizonless objects and their stability properties. As explained in the Intro¬ 
duction, these objects are implied by unitarity in lien of black holes and, for 
the sake of brevity, are also referred to as stars here. They are not stars in a 
conventional sense because, for example, their central densities are likely to 
be of the order of Planckian densities; and their sizes are likely to be of the 
order of Schwarzschild radii. 

The spacetime is assumed to be D = n c + rn + 2 dimensional with rn > 2 , 
with n c dimensional compact toroidal space, and with (rn + 1) dimensional 
non compact space. The stars are assumed to be static and spherically 
symmetric in the non compact space, and to be made up of non interacting 
multicomponent fluids with linear equations of state. The M theory stars 
are eleven dimensional and are assumed to be made up of A f stacks of M2 
and M 5 branes, intersecting according to the BPS rules whereby two stacks 
of five branes intersect along three common spatial directions; a stack each 
of two branes and five branes intersect along one common spatial direction; 
and two stacks of two branes intersect along zero common spatial direction. 
Assuming the spatial directions of the brane worldvolumes to be toroidal 
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and assuming necessary isometries, the intersecting M theory branes can 
be described by A f seperately conserved energy momentum tensors T M n(i), 
I = 1,2, • • •, J\f , with appropriate equations of state among their components 
[55, 59, 60, 61, 56, 57, 58, 63, 64, 65], 

Following closely the notations and conventions of our earlier work [37], 
we now write a suitable ansatz for the metric and obtain the equations of 
motion. Let x M = (t, x 1 , r, 6 a ) be the spacetime coordinates where x l , i = 
1,2, • • • ,n c , describe the n c dimensional toroidal space; and the radial and 
the spherical coordinates (r, 9 a ), a = 1,2, describe the (m + 1) 
dimensional non compact space. In standard notation and with n 2 = 8 tiGd = 
1 , the equations of motion may be written as 

TZmn ~ ~ 9mn 77 = Tmn = T Tmn(i) (1) 

2 i 

£ v « T %> = 0 (2) 
M 

where Tmn is the total energy momentum tensor of non interacting mul¬ 
ticomponent fluids and Tmn(I) is the energy momentum tensor for the I th 
component fluid. For stars in M theory, Tmn is the total energy momentum 
tensor for intersecting branes and Tmn(i) is the energy momentum tensor for 
the I th stack of branes. 

In the following, we consider static solutions which are spherically sym¬ 
metric in the (m + 1) dimensional non compact space. We will study the 
singular solutions and the asymptotic perturbations around them in the limit 
of large r [34, 29]. The suitable ansatz for the line element ds is given by 

ds 2 = gMN dx M dx A ' = —e 2X °dt 2 + ^ e 2y {dx 1 ) 2 + e 2X dr 2 + e 2<T dh2^ (3) 

i 

where dVt m is the standard line element on an m dimensional unit sphere. 
The energy momentum tensors T M n(i) are assumed to be diagonal. These 
diagonal elements are denoted as 

(n { 7 ), T r r (I)i T a a (I)) = (PON Pill n 7 , p aI ) 

where p 0 i = — Pi and p a j = pr for all a . The total energy momentum tensor 
is now given by T M N = diag ( po , pi, II, p a ) where po = — p , p a — p for all 
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a , and 


p = J2pi > Pi = Ep^ > n = ^n 7 , p = J2pi 


Define 


a = (0, z, a) , A" = (A, A*, A a ) , p aI = (p 0 i , pu,p a i) 


where A a = cr for all a . Also, define 

a = E a “ = a ° + E ai + ^ > T I = '£T% I) = n I + 'EPai • 

a i M ol 

For static solutions which are spherically symmetric in the (m + 1) dimen¬ 
sional non compact space, the fields (A", A) and ( Poi i, If/) depend only on the 
coordinate r . Using the above definitions and the metric given in equation 
(3), it follows straightforwardly that the equations of motion (1) and (2) now 
give 


(n,)r 

= -n 7 A r + J2p«iK 

Ol 

(4) 

£(V) 2 

= 2 Y, Hz e 2A + m(m - 1) e 2A " 2CT 

(5) 


a I 


K r + (Ar - A r ) A“ = J2 (-p Q j + e 2A + (m - 1) e 2A “ 2CT (6) 

where the subscripts r denote r— derivatives. We also define a function /(r) 
and a mass function M(r) by 


„2A-2cr 


1 

r 2f 


/=!- 


M 

~m— 1 


so that either of them may be traded for the function A(r) 


Reduction to d = m + 2 dimensions 


We will now dimensionally reduce on the n c dimensional toroidal space 
from D dimensions to d — m + 2 dimensions described by the x M = (t, r, 9 a ) 
coordinates. Consider the D dimensional line element ds given by equation 
(3), and denote its d dimensional part as follows: 

ds 2 d = gijjy(d) dx^ dx u = — e 2X °dt 2 + e 2X dr 2 + e 2a dVL 2 n . 
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Upon dimensional reduction, symbolically, we have 

S ~ j d D x R ~ J d d x yj-g(d) e A ° (R (d ) H-) 

~ J d d x (R H— •) 

2A C 

where A c = X; A* and = e~ g^{ d ) is the d dimensional Einstein frame 
metric. The corresponding line element ds d then becomes 

ds 2 d = g, w dx^ dx v = —e 2X °di 2 + e 2X dr 2 + e 2a d£l 2 m (7) 

where 

A c - A c 

A Q = A Q H-, A = A H-, (8) 

m m 

and A a = a for all a . Furthermore, one has 

A = A 0 + md = A + — , e 2A " 2 " = e 2A - 2ff = • 

m r z j 

Note that A* and A* can be expressed easily in terms of each other 8 and, 
hence, they are both equally convenient to work with. 

Writing p 0I = —pi, it now follows from equations (4) - (6) that 

( q- \ 2U 

Pa -- ) A* H-- A c r (9) 

m / m 


2A ° r a r + (m — l)(oy) 2 


cr rr | ^A..p Aj.^ <7 r 

X° rr + (A r - A r ) X° r 
K r + (Ar - X r ) At 


-£ll 7 e 2A + (m — 1) e 2A - 2 " + - 


rn 


T r 


~Pr H-e 2A + (m - 1) e 

' m ' 


m 

2A-25- 


E(ft + - 

j V m 

{-'Pit + 


„2A 


Tr 


m 


2A 


( 10 ) 

(11) 

( 12 ) 

(13) 


8 For any a a s, let a a = a a + ^ where a c = X» a * • Then a c = X» = ( n c + m )m and, 
hence, a a = a a — n a +m . Similarly for any b a s. Also, X* = X» a 4 6* — X» 
which further implies that X* > 0 if a 1 do not all vanish. 
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where Ti — Hi — pi + m pi and B = J2i K K ■ Using the diffeomorphic 
freedom in defining the radial coordinate, we now set e a = r . Equations 
(10) and (11) become 


r X° r 
r (A° - K) 


r X r 


E 

i 

E 

i 

E 


E , ,, m — 1 

— r 2 e 2A H- 

m 2 


(e 2A - 1) + 


r 2 B 


fU - pi 


m 


Pi „2 „2\ m 1 /„2\ 


2m 

r 2 e 2A + (m - 1) (e 2A - 1) 
r 2 B 


re - 

m 2 


(e 2A - 1) + 


2 m 


(14) 

(15) 

(16) 


Since e a = r , we also have 


ds 2 d = 


^°de + ^-+r 2 dSl 2 m 


f = 


-2A 


1 - 


M 

rplTL— 1 


Thus, the line element ds d in the d = m + 2 dimensional Einstein frame takes 
the standard form. Therefore the functions appearing in it, e.g. the mass 
function M, may be interpreted in the standard way. For instance, Madm, 
the ADM mass of the star of radius r* is related to the mass function by 

ra+1 

M (r*) = 16 J r (t ,P Madm where S m = 2 J m +i . is the ‘area’ of an rn dimensional 

x 7 BIX Dm Vn c I (—^— ) 

unit sphere and V n<: is the coordinate volume of the n c dimensional toroidal 
space. Note that equation (16) and the relation M{r) = r m_1 (1 — e _2A ) 
now give 

Mr = E- . (17) 


Linear equations of state 

To solve equations (9) and (12) - (15), and to obtain solutions for the 
fields (A Q , p a i , II/), one further requires equations of state which give p a j and 
11/ as functions of p T . In cosmological contexts, it is a standard practice to 
use linear equations state. In the contexts of stars, Chavanis explains several 
situations where such equations of state arise and the reasons for using them 
[29]. In the case of horizonless objects which are implied by unitarity and 
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which are of main interest here, nothing is known rigorously about the nature 
and properties of their constituents. Hence, it is not presently possible to 
derive the equations of state from the underlying microscopic physics. In 
order to make progress, we will assume here linear equations of state. The 
stiffest equation of state (pressure = density, sound speed = light speed) 
is of this type. Also, the linearity helps in obtaining explicitly the singular 
solutions and the asymptotic perturbations around them, and in studying 
the stability properties. 

Thus, let 

Pal = W I a Pl , bl/ = w[ pi (18) 

where w J a and ud are constants, w : 0 = — 1 since p 0 j = —pi, ud = ud since 
Pal = Pi f° r all a , and we assume that ud > 0 . It is common to take 
the pressures inside the stars to be isotropic along the noncompact spatial 
directions, namely to take the pressure H/ along the radial direction to be 
equal to the pressure pr along the transverse spherical directions. However, 
the assumption about isotropy may be unwarranted and restrictive : it is not 
required by spherical symmetry and, in general, the pressures inside the stars 
may be different along the radial and the transverse spherical directions of 
the non compact space. For conventional boson stars, anisotropic pressures 
arise naturally when scalar fields are present [38, 39]. See also the reviews 
[40, 41, 42, 43] and the recent paper [44], We will assume here that such 
an anisotropy may be present and, hence, that ^ pj in general. As a 
measure of this anisotropy, we define a dimensionless parameter rj 1 by 

r m fHi — pi\ m f ud — w T ^ 

11 = ~2 V n 7 ) = ~2 V K ) 

so that rj r = 0 corresponds to the isotropic case. The factor of y is for 
convenience. 

Since Ti = Hi — pi + m pi, equations (18) give 

(~PaI + -)= C aI Pi , C aI = +W I + . 

V m / m 

Hence 

c 0/ = 1 + w 1 + c 1 , c u = — «,d + W 1 + c 1 , c aI = c I = . (20) 

m 
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, , V. e*' 

The corresponding untilded coefficients are given by c = c — n J +m , see 
footnote 8. Let 0 7 be given by 

= -(l + ^)A°-2 f ,ya-^c t/ A I . (21) 


It then follows from equations (9) and (19), and from e a = that 


Pi = pi o e*' e“ , r-pi e 2 "' = p; 0 e*' +2A+2 ‘' 


( 22 ) 


where p/ 0 > 0 is a constant. One also obtains, for any function A(r(A)), 


Xpj = rX r , A'a-g- = r 2 X rr + rA r 
?’ 2 + (A. r — A r ) X^j = Xm + (xa- — A^) X& 

where the subscripts a denote a —derivatives and x — A — a = A° + (m — 1) cr . 
Equations (12) - (16), written in terms of A, now become 

AL + (X* - A.) A° = E 7 07 PI0 e <t >I+2 *+ 2a (23) 

i 

AA + (x* ~ A*) Al = E A 7 PI o e^ +2 * +2 * (24) 


A 


A? - A, 


A, 


E — P/o e^ +2 * +2 * + 

j m 2 

E 5 7 P/o e^ +2 * + (m - 


V- AfO <p+2A+2<r m 1 

Y m 2 




| r 2 ff 
2m 


(25) 

(26) 

(27) 


where ?’ 2 15 = A| A-. . The above equations thus describe stars whose 
constituents obey the linear equations of state (18). 


4. Asymptotic analysis : Singular solutions and perturbations 
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Consider the limit r = e a —> oo . Suitable ansatzes for A 0 , A*, and A in 
this limit are given by 

A 0 = s° a + u° , A 1 = s' 1 a + u l , A = A 0 + s a + u (28) 

where s°, s\ s, and A 0 are constants and u°, u\ and u are functions of r . 
Further constants A[] and A ? 0 could have been added to A 0 and A* also but, 
with no loss of generality, they have been set to zero. Writing 

(f) 1 = q 1 o + y 1 (29) 

it follows from equation (21) that q 1 and y 1 are given by 

wi q 1 = - (1 + «/) s° - 2 r/wl - ^ c u s* (30) 

i 

wW = . ( 31 ) 

i 

Also, write r 2 B = E; K At = B 0 + 2Bi + B 2 where 

e„ = ys i s ‘, b 2 = E44- 

iii 

In the limit r —>■ oo , the Ao and the a terms in the above equations give the 
leading zeroth order asymptotic solutions. They give the singular solutions 
of [34, 29]. The functions u 's are treated as perturbations and are used to 
obtain the first order corrections to the leading asymptotic solutions. They 
give the perturbations around the singular solutions. 

Zeroth order : Singular solutions 

Consider the equations of motion (23) - (27) and expand them to zeroth 
and first order in the functions it’s. At zeroth order, equating the powers of 
r gives 

2 + q 1 + 2s = s = 0 (32) 
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2 . 9 Now, equation (30) becomes 


and, hence, q 1 = — 


2 wi (1 - ri 1 ) = (1 + wi) §° + £ c u s 1 . ( 33 ) 

i 

Also, upto first order in the functions u’s, we have 

r 2 ( Pl e 2X ) = p 70 e 2W+2 “ = R I (l + y I + 2u+ •••) 

e 2A _ i = e 2A 0 +2 u _ x = ( e 2A 0 _ 1) + g2Ao ( 2 U+ ■■■) 



a + («~ — ug) , a = in — l + s° — s 


where Rj 


Pro e 2X ° . At zeroth order, equations (23) - (27) now give 


as = 


= E'c 0I Ri 


~s° = 


a s 


= £ i“ R, 

I 


= J2 — R i + 


rn — 1 


(e 2 ^ 0 - l) 


+ 


Bo_ 

2m 


s — s = 


Y^c 1 Ri + ( m — 1) ( e 2X ° — l) 


s = 


= 


Ri m — 1 


rn 


(e 2Xo - l) 


1 ' + ^- 


(34) 

(35) 

(36) 

(37) 


and equation (36) gives 

a = m — 1 + s° — s = ^2 c 7 Ri + (m — 1) e 2 ^ 0 . (38) 

i 

Equations (33) - (38) constitute the equations of motion at the leading zeroth 
order in the functions u’s. They will give the singular solutions. 


First order : perturbations around singular solutions 

9 In general, one should analyse equation (30) which determines q 1 and thus the asymp¬ 
totic behaviour of pi . For a given set of values for and w 1 , some of the resulting q 1 s 
may lead to subdominant terms. Then the corresponding pi s become unimportant and 
effectively reduce AT . With no loss of generality, we are assuming that ud and w 1 are 
such that q 1 = —2 for all /, thus all pi s remain important and AT remains unreduced. 
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At first order in the ft’s, equations (23) - (27) give 


uL + au l + s° 0° - ua) = Y 507 R i ( y 1 + 2u ) 


Ka + a K + s* {u% - u 5 


= Y ? 1 Ritf + Zu) 


u 


° = £ — Ri {y 1 + 2 u) + (m — 1 ) e 2X ° u+^ 1 

j m 


m 


w~ — Us- = £ c 7 R i ( V 1 + 2 m) + (m — 1 ) e 2X ° ( 2 u) 

1 

= £ — {y 1 + 2 u) — (m — 1 ) e 2 ^ 0 u + — . 

j m m 

Using the zeroth order results for the u —terms in the right hand sides of the 
above equations, one obtains 


m~ = Y — R i y 1 + (2s° + m — 1 — —u + — (39) 


m 


m 


m 


2 °. 

l aa 


-u- — = £ c 7 Rj y 1 + 2cm 

l 


(40) 

sr' R j 1 ( i 

Ua = Y — y - I m- H- 

j m V m) 

u-\ -. 

m 

(41) 

and the u l ad equations now become 



n^ + aul = £(c 07 -m 0 c 7 ; 

7 

) Ri y 1 

(42) 

u'tt + au'z = £(£-§£) 

Ri y 1 

(43) 


Equations (31) and (39) - (43) constitute the equations of motion at first 
order in the functions u’s. Their solutions give the perturbations around the 
singular solutions. 


Significance of the perturbations 
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We now outline briefly the significance of the singular solutions and the 
perturbations around them. In section 2, we have briefly mentioned a few of 
these aspects. See [29] for complete details. 

The singular solutions and the perturbations around them can be used, 
among other things, as indicators of the stability of stars. A star, with its 
constituents obeying linear equations of state, is enclosed in a box of radius 
r* so as to prevent the evaporation of its constituents and to make its total 
mass finite. The mass of the star is then given, upto constant numerical 
factors, by the mass function M (r*) evaluated at r* . The singular solutions 
give the scaling relations between mass, radius, and other quantities of the 
stars. The singular solutions are of no relevance when the stars are unstable 
but, when stable, these solutions give the scaling relations in the limit of 
large central density or large radius of the stable stars. 

The behaviour, namely damped oscillatory or monotonic non oscillatory, 
of the perturbations around the singular solutions leads to the corresponding 
behaviour for the mass - central density profile in the asymptotic regions. 
Considering a series of equilibria, Chavanis obtains the mass - central density 
profile in the asymptotic and non asymptotic regions. He shows that equi¬ 
librium configurations become unstable beyond the first maximum in this 
profile. This is also shown to correspond to conditions for nonlinear dynamic 
stability. In those higher dimensional cases where the behaviour of the per¬ 
turbations is monotonic non oscillatory, Chavanis obtains the mass - central 
density profile in the non asymptotic regions also and shows that the entire 
profile is monotonically increasing and asymptotes to the constant value line 
given by the singular solutions. 

In this paper, assuming linear equations of state, we will obtain singular 
solutions and perturbations around them. Enclosing the star in a box of ra¬ 
dius r*, the mass - central density profile and its behaviour in the asymptotic 
regions can also be obtained from the perturbations. However, in the case 
of horizonless objects of interest here, nothing is known rigorously about the 
nature and properties of their constituents. Hence, we are unable to carry out 
the analogs of the analyses mentioned in the previous paragraph. Namely, 
since the detailed properties of the constituents are not known, we are unable 
to obtain the profile in the non asymptotic regions, and to study the stability 
properties around the maxima. 

We note below a few useful points. 
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(1) The expressions = 1 — e 2A and A = Ao + u imply that, in the 
limit of large r, 

= 1 - e _2A ° + 2 e _2A ° u + ■ ■ ■ . (44) 

y>m—l v ' 

The first two terms correspond to the singular solutions and the u term 
describes the perturbations in the mass function. Similarly, y 1 describes the 
perturbations in the density p 1 . 

(2) The radius r* of a spherically symmetric star, defined to be given 
by II(r*) = 0, is infinite when it is made up of perfect fluids with linear 
equations of state, see [33] for a derivation. Hence, let the star be enclosed 
in a box of radius r* which will render its mass 

M(r*) ~ J dr M r 

finite where M r is given in equation (17). Then, following the analysis of 
[29], the perturbations around the singular solutions can be used to obtain 
the mass - central density profile of the star in the asymptotic limit of large 
central density or large radius. 

Let x c h oc y rjTc r* be a measure of central density p c , and let y c h = 
be a measure of the mass of the star. By detailed analysis of the equations, 
and incorporating the properties of the constituents of the star, Chavanis 
obtains the mass - central density profile in both the asymptotic and non 
asymptotic regions and finds that : (i) As x c h increases from zero to oo, 
y c h increases from zero to a (first) maximum y\ at aq, thereafter exhibits 
damped oscillations, asymptoting to a value y s . (ii) The behaviour for large 
values of x c h can be seen from the singular solutions and the asymptotic 
perturbations around them, (iii) The solutions are unstable beyond the first 
maximum which is at ( x c h , y c h) = (aq, yi) . 

As a consequence, one has the following. For a given value of central 
density, the radius r* must be < rq* where rq* corresponds to aq. The mass 
of the star must then be less than (j/i r™ _1 ) . A more massive star will be 
unstable and will collapse. 

(3) For D = m + 2 dimensional stars, n c = 0 and J\f = 1 in our nota¬ 
tion, Chavanis finds that if m > m cr ~ 9 then y c h increases monotonically 
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from zero to y s , effectively making x\ infinite and y\ — y s . See Figure 
23, and also Figures 20 - 22, in [29]. The asymptotic perturbations around 
the corresponding singular solutions exhibit monotonic behaviour with no 
oscillations. 

As a consequence, one has the following. For m > 9, X\ is effectively 
infinite. Then, for a given value of central density, rq* is infinite which makes 
the upper limit (y s r™ _1 ) on the mass of the star also infinite. Hence, a star 
can be arbitrarily massive and stable when m > 9 . Also, in the limit of large 
central density or large radius, one obtains the scaling relation Af(r*) ~ r™ _1 
which shows that the sizes of the large stable stars scale as their Schwarzschild 
radii. 

(4) In [37], we generalised this study to D = n c + m + 2 dimensional 
stars, with n c toroidal directions, made up of AT > 1 number of perfect fluids 
where pressures are isotropic along the non compact spatial directions. Stars 
in M theory correspond to specific values of n c and A f . We found that, even 
in these generalised cases, m > 9 is required for stability. 

5. Singular solutions and asymptotic perturbations 

Taking the pressures inside the stars to be anisotropic along the radial 
and the spherical directions of the non compact space, and hence taking 
the anisotropy parameters y 1 ^ 0, we now obtain singular solutions to the 
equations of motion and asymptotic perturbations around them, generalising 
those given in [34, 29, 37]. These solutions may be obtained for any general 
set of values for n c , m, A7, and w^. However, such a generality is 

neither illuminating nor needed for our purposes here. Hence we will present 
only three cases which are illustrative and are also of direct interest. 

In this section we will consider (in + 2) dimensional stars made up of a 
single fluid. In the next section we will consider two examples of stars in 
M theory : In one, the stars are made up of a stack of M2 or M 5 branes. 
In another, the stars are made up of four stacks of intersecting branes, two 
stacks each of M2 and M 5 branes, see footnote 5 . 

(m + 2) dimensional stars with J\T — 1 

Consider (m + 2) dimensional stars made up of a single fluid. Then 
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n c = 0 and AT = 1 . The tildes and the /—scripts on various quantities are 
now unnecessary. Hence we omit them. 

The fluid is assumed to have anisotropic pressures in general. Its radial 
and the transverse spherical pressures, and the anisotropy parameter q , are 
given by 

m fw n — w 

H = w n p , p = w p , q = — - 

2 V w n 

Also, c° = 1 T w T c and c a = c = Wn ~ 1 and hence 

m c° = (m — 1) (1 + w n ) + 2 w n (1 — q) 



We now write down the leading order asymptotic solutions to the equations 
of motion and perturbations around them. We have s = 0 and q = —2 . The 
zeroth order equations (33) - (38) then give the following relations. 


R 


(in — 1 ) e 2A ° 


2 uy (1 — g) 

1 + uq 

a s° 2 m Wn (1 — q) 

C° (1 + Wn) 2 

V 

a — c R = -- 

(l + O 2 


a = m — 1 + s u 


m c° 

1 + Wn 


V = (m- 1) (1 + Wn ) 2 + 4 uv (1 - q) 


The above expressions describe the singular solutions, 
give 


e 2A ° ~ r 2s ° 


2\ 




AiUn(l~q) 

( 1+Wn ) 2 


For example, they 


Po 2 m(m — l)wn(l ~ q) 

P — r 2 — r 2 V 

Note that, from R oc p 0 > 0 and w n > 0 , it follows that q < 1 and, hence, 
that me 0 , s°, a, and V are all positive. The mass function can be obtained 
from equation (44) and, using the above expression for e 2A ° , it is given by 


^T = § ( 2 w * ( 1 ~V) + (' m -!)(! + Wir) 2 « + •••) • (45) 
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The mass of the star M(r*) enclosed in a box of large radius r* is then given 
by 

Mir.) * C- • 

Consider now the first order equations of motion (31) and (39) - (42). 
Their solutions will give the asymptotic perturbations around the singular 
solutions. Equation (31) gives, upon using the expressions for s° and R , 

u° 

s° V = ~ 2 (1 - 77 ) u° , Ry = - 2 (1 -r)) a — . 

Equations (39) and (41) give 

u° a = —s° u° + (m— 1 + 2s°) u 

s° 

u a = - u° — (m — 1) u 

w n 

from which it follows, after a little algebra, that both u° and u obey the same 
equation given by 

(»U + a (*)„ + (t) = ° (46) 

where (*) = u° or u . ft is straightforward to show that equation (42) also 
gives the above equation for u° . 

The solutions to equation (46) are of the form (*) ~ e k<7 where 

k = ~ a±x/K . A = a 2 _ 8 ( 1 -’?)P . 

2 (l + O 2 

If A < 0 then, in an obvious notation, k — — k re ± i k im with k re > 0 . 
The solutions for (*) are then oscillatory and, since a = In r , they may be 
written as 

A 

(*) = ~~j~~ Sin (k im In r + B *) 

rpfcre 

where A* and B * are integration constants. If A > 0 then a > \fK since 
1 — 77 > 0 and V > 0 and, in an obvious notation, k = —k\ ± k 2 with 
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h> k 2 >0 . Hence the solutions for (*) are non oscillatory and they may 
be written as 

If A = 0 then k 2 = 0 and the solutions are (*) = 4y (In r + B*) and are non 
oscillatory. Considering the star enclosed in a box of radius r*, introducing 
Milne variables, and following the analysis of Chavanis given in [29], one can 
now obtain the mass - central density profile in the asymptotic limit of large 
central density or large radius. The asymptotic behaviour of this profile is 
similar to that of the perturbations: it is oscillatory if A < 0, and is non 
oscillatory if A > 0 . 

Comparing with the work of Chavanis in [29], we note that the solutions 
for (*) given above when A < 0 and when A > 0 will lead to the analogs 
of equations (174) and (175) in [29]. The resulting mass - central density 
profile for M(r*) will lead to the analogs of the asymptotic parts of Figure 
23 in [29]. The numerical and the analytical studies leading to the analogs 
of the non asymptotic initial parts of that Figure are beyond the scope of 
the present paper and hence, although important, are not attempted here. If 
one assumes that these non asymptotic initial parts remain qualitatively the 
same for the anisotropic case also, then one may conclude that a star can be 
arbitrarily massive and stable when A > 0 ; and that the singular solutions 
describe its mass - radius relations in the limit of large central density or 
large radius. 

under which A < 0 and A > 0 . Let 


A = A wl + 2B w n + C (47) 

where, after some algebra, one obtains 

A = (m — 3 + 2 if ) 2 

B = (m + 1 — 2rj) (m — 9 + 877) 

C = (m — 1) (m — 9 + 877) , 


We now study the conditions 


A = 


{X-^W'K ) 


2 . It then follows that 
—m c° ± \fK 

k = 


2 (1 + w n ) 
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AC — B 2 = 32 (1 — rj) 2 (m — rj) (m — 9 + 8 rj) . 

Consider now the sign of A . It is same as that of A . Using the quadratic 
expression for A given above, it can be seen that if AC — B 2 < 0 then A can 
be negative for a range of values for m, rj, and iv n . For example, if m = 2 
and 77 = 0 then AC — B 2 < 0 . Then A = w 2 — 42 w n — 7 and is negative, for 
example, for 0 < w n < 1 . If AC — B 2 > 0 then A > 0 always. 


In the isotropic case, rj — 0 and the above expressions reduce to those 
given in [29]. Then A > 0 if m > 9 and, depending on the value of w n = w , 
A(w) > 0, if m > m cr {w ) ~ 9 . In the anisotropic case, rj is non vanishing. 
Noting that ij < 1 and m > 2 , it follows from the above expressions that 
AC — B 2 > 0 if the values of rj and lie in the range given by 


9 — m 
8 


< rj < 1 


5 m — 9 w 
Am ~ w n 


m — 2 
m 


(48) 


then A > 0, hence A > 0, and the solutions are non oscillatory. The de¬ 
pendence of these ranges on uq can also be incorporated, but is superfluous 
for our present purposes of showing that anisotropy can lead to asymptotic 
non oscillatory solutions. Note that the isotropic case rj — 0, equivalently 
w = w n , is included in the ranges given above only when m > 9 . For 
lower values of m, (for example, for m — 2 which corresponds to four di¬ 
mensional spacetime) certain amount of anisotropy ( rj > | ) is needed to 
obtain the non oscillatory behaviour of the asymptotic perturbations around 
the singular solutions. 


6. Stars in M theory 

In this section we will analyse stars in M theory. In section 2, we have 
described briefly M theory stars and how the necessary equations of state are 
obtained. We now proceed with the analysis. 

The spacetime is eleven dimensional in M theory, having n c dimensional 
compact toroidal space and (m + 2) dimensional non compact spacetime 
where n c + rn = 9 . The M theory stars are taken to be made up of A f 
stacks of M2 and M 5 branes, intersecting according to the BPS rules. We 
take the spatial directions of the brane worldvolumes to be toroidal and 
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assume necessary isometries. These intersecting branes can be modelled by 
J\f number of seperately conserved energy momentum tensors [55, 59, 60, 61]. 
Hence, the present formalism can be applied to the corresponding stars. 

M theory has U duality symmetries. As shown in [58], see also footnote 
7 , they lead to a relation among the components (p a i, Hj) of the energy 
momentum tensor Tmn(i) which is given by 

P\\i = n / + p 0I + p±i + m (p 7 - p_ u ) 

where pp and p±i are the pressures along the directions that are parallel and 
transverse to the worldvolume of the I th stack of branes. The above relation 
is a consequence of U duality symmetries and, therefore, must always be valid 
independent of the details of the equations of state. Also, since the sphere 
directions are transverse to the branes, it is natural to set pj = p±j . 

Note that, for stars in M theory, there is no compelling reason to take the 
pressures H/ and pj to be equal. Indeed, in the case of charged intersecting 
black branes in M theory, these pressures are not equal although the relation 
Pi = p±i and the U duality relation above are obeyed. Hence we assume 
that, in general, H/ ^ pi for stars in M theory. Setting p±i = pi and 
Poi = —pi , the U dulaity relation now becomes 

P\\I = - Pj + pj . (49) 

Consider the linear equations of state given by (18). Let pp = wf pi and 
P±i = w{ pi where wf_ = w 1 since p±i = pi . The LI duality relation then 
gives 

wf = wf — 1 + w 1 . 

Now consider the coefficients c rl = —wf + W 1 + C 1 defined in equation (20). 
Note that wf — if i e ||/ , namely if x 1 is a worldvolume coordinate of the 
I th stack of branes; otherwise, wf = wf_ . It then follows from iuf_ = iv 1 and 
the LI duality relation for wf given above, that the corresponding coefficients 
c" 7 and cr 17 are given by 

c 117 = — wl + 1 + c 7 = (1 — m) c 1 , c ±7 = c 7 = —-- . (50) 

m 

Hence, for any a 1 with a 1 = a 1 + fff and a c = a* , it follows that 

(51) 

i i i&\\i 
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Y, c 11 a 1 = Y ~ ciI ai = c 1 a c - m c 1 Y a * 



In passing, we note that the relation c^ 1 = (1 — m) c^ 1 is same as that 
obtained in the corresponding isotropic cases studied in [37]. Consequently, 
the resulting equations for M theory stars will have very similar structure in 
both the isotropic and the anisotropic cases. With ud and c aI specified for 
intersecting branes, we now consider two examples of stars in M theory. 


M theory stars made up of M2 or M5 branes 

Consider stars in M theory made up of a stack of M2 or M 5 branes. Then 
AT = 1, n c — 2 or 5, and m = 7 or 4 . The /—scripts on various quantities 
are now unnecessary and, hence, we omit them. Also, we will first write the 
solutions in a form applicable for J\f = 1 and for any values of n c , m, and 
c l ; and then, at the end, specialise to the case of M2 or M 5 brane stars. 

We now write down the leading order asymptotic solutions to the equa¬ 
tions of motion and perturbations around them. We have s = 0 and q = —2 . 
The zeroth order equations (33) - (36) then give the following relations. They 
describe the analogs of singular solutions in this context. 


2 W n (1 — 7}) = (1 + W 7 ,) S° + £ C * S l 

i 

-° - c° R , Q? = c* R 


as = 


s = 


s = 




m— 1 


R + 9 

m 2 


(e 2 ~ Xo - l) + 


B [o_ 

2m 


c R + (m — 1) (e 2 ^ 0 — l) 
Using these relations, one obtains 


R = a 


s l = R 


~ 0 = 2 w* (1 - V) 


(1 + w n ) (1 + 7) 


y c 1 e 

where 7 = y - 0 ; and, after some algebra, 


m — 1 ) e 2A ° = a — c R = 


a 


V 


1 + 7 c° (1 + uq) 


+ 7 


(52) 


(53) 
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where V — (m — 1) (1 + w n ) 2 + 4% (1- rj) . Note that, from R oc p 0 > 0 
and w n > 0 , it follows that rj < 1 ; and, hence, that me 0 , 7 , s°, a, and 
T> are all positive. Note also that the effects of the compact toroidal space 
appear through the parameter 7 alone. 

Consider the first order equations of motion (31) and (39) - (43). Their 
solutions will give the asymptotic perturbations around the singular solu¬ 
tions. Upon using s l = c l % , equations (42) and (43) give 

F-~ + a F- = 0 

where F l — u l — c* 4^ . Although it follows that F l = F[ e~ aa in general, 
we will set the integration constants F\ to zero. This gives u l = c l , which 
we now use in the remaining equations (31) and (39) - (42). 

Equation (31) gives, upon using equations (52) for 5° and R , 

u° 

S° y = - 2 (1 - 77 ) u° , Ry = - 2 (1 - 77 ) a — . 

Equation (39) gives, after some manipulations, 10 

u° a = - (1 + 7 ) 5° u° + (m - 1 + (2 + 7 ) 5°) u . 

Equation (40) gives straightforwardly 

u° 

ul — U& = — 2(1 — rj) ac — + 2 a u . 

c u 

After a little algebra, it follows from these two equations that both u° and u 
obey the same equation which we write as 

(*)?? + « (*)* + 2 (1 -rj)(a-c R) (*) = 0 (54) 

where (*) = u° or u and (a — c R) is given in equation (53). Equation (42) 
gives the equation for u° straightforwardly in the above form. 

The solutions to equation (54) are of the form (*) ~ e ka where 

k = — ° ^ , A = a 2 — 8 (1 — 7 ) (a — c R) . 

10 We used the expressions B 0 = 7(1 + w v )^-, Bi = 7(1 + w„) u% , and = 
a + 7 s° , which can all be derived easily. 
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As explained in the case of (rn + 2 ) dimensional star, it is important to study 
the sign of A since it determines whether the solutions for (*) are oscillatory 
or not. This behaviour of the perturbations (*) will lead to corresponding 
asymptotic behaviour in the mass - central density profile of a star enclosed 
in a box. Writing 

A = a(A + £) 

m 5° (1 + w w ) (l + y) 

it can be shown after a long but straightforward algebra that A and 8 are 
given by equation (47) and 

8 = m (m — 9 + 87 ) ^ c l c l . (55) 


Consider the sign of A . It is same as that of (A + 8) . If c 1 = 0 for all 
i then 8 — 7 = 0 and the earlier analysis given below equation (47) applies 
directly. Now let c* do not all vanish. Then cV > 0 , see footnote 8 . It 
then follows that 8 > 0 if m — 9 + 877 > 0 . Hence, noting that 7 < 1, one 
has 8 > 0 if the inequalities (48) are satisfied. Interestingly, this is precisely 
the condition that also ensures that A > 0 but we do not know a simple 
reason, if any, for this coincidence. We thus have that A > 0 and the 
solutions are non oscillatory if the values of the anisotropy parameter 7 and, 
equivalently, of — he in the ranges given in equation (48). Also, note that 
the value 7 = 0 , equivalently w = w n , is included in these ranges only when 
m > 9 . Otherwise certain amount of anisotropy is needed along the non 
compact spatial dimensions to obtain the non oscillatory behaviour of the 
perturbations around the leading order asymptotic solutions. 

The solutions obtained above are for J\f = 1 and are applicable for any 
values of n c , m, and c l . In particular, they are also applicable to stars in M 
theory made up of a stack of M2 or M 5 branes for which n c = 2 or 5, m = 7 
or 4 , and c* = = (1 — m) c as given in equation (50). It then follows easily 


that 


£cv 


n c m 
n c + m 



A nc {w n - l ) 2 


where A Uc = | for M2 branes and = | for M 5 branes. 


M theory stars made up of J\f = 4 , 22'55' intersecting branes 
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Consider stars in M theory made np of two stacks each of M2 and M 5 
branes intersecting according to the BPS rules whereby two stacks of five 
branes intersect along three common spatial directions; a stack each of two 
branes and five branes intersect along one common spatial direction; and 
two stacks of two branes intersect along zero common spatial direction. Let 
22'55' : (12, 34,13567, 24567) denote the configuration of the intersecting 
branes and indicate the spatial worldvolume directions of the four stacks of 
branes. Now Af = 4, n c = 7, and rri — 2 . Also, the coefficients c l1 are given 
by equations (50). Thus, 



Note that equation (51) gives, with S c = J2i an d U c = u% , 


£ c a s' 

= £ c u s‘ 

= c 1 S c — m c 1 ^2 s l 

(56) 

i 

i 

tell / 



- £?'«• 

= c 1 U c — m c 1 ^2 u z . 

(57) 


i i »e|U 


We now write down the leading order asymptotic solutions to the equa¬ 
tions of motion and perturbations around them. They will describe the 
analogs of the singular solutions and the asymptotic perturbations around 
them. We have s = 0 and q 1 = —2 . Equation (33) can then be satis¬ 
fied for all I by choosing iu^ = w w and w 1 = w for all I , see footnote 9 . 
Then q 1 = q, c 1 = c, and c oi = c° for all /, but c lT do depend on / since 
c l1 = if i E ||/ and c l1 = c ± otherwise. Now, using equation (57) for the 
22 , 55 / : (12,34,13567,24567) configuration, note that 

(Z cil = c U c {J\f - 2 m) = 0 (58) 

since Af = 4 and rri — 2 . 11 Consider sums of the type XA c l/ Xj . Upon 

11 Cancellations of this type do not happen for all intersecting brane configurations. It 
happens in the present case, and for the case where three stacks of two branes intersect 
(for which Af = 3, n c = 6, and m = 3), and for the equivalent U dual versions of these 
two configurations. Similar cancellations happen for these two independent configurations 
in the cosmological context also. There, the cancellations may be understood as arising 
due to the balancing of contraction or expansion forces applied by the branes on the 
compcat directions parallel or transverse to the worlvolume directions [63, 64, 65]. These 
cancellations are responsible for the stabilisation of the compact toroidal directions. 
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using fill = — c 1 - and denoting the X/s as (X 2 , XV, X 5 , X 5 /) , such sums 
become 


E 


E^Ab 

T 

— gll 

E 

to 

1 

* 

to 

+ 

* 

Cn 

1 

* 

oj 

"M 1 

to 

* 

>—1 

= c 11 

(X 2 - X v - X 5 + X 5 0 

E5 37 Ab 

/ 

= c 11 

(~X 2 + Xy + X 5 - x 5 ') 

E^Ar 

/ 

= c 11 

(-X 2 +Xy-X 5 + Xy) 

C 5 ’ 6 ’ 7 7 X] 

= c 11 

(-X 2 -X 2 , + X 5 + X 5 ,) 


Consider equation (33). Upon using equation (56), it gives 

V- J (1 + w n ) S° - 2 w n (1 -rj) i S' 0 

/ ^ & ~ I 1 

zrf me m 


(59) 


which implies that the sum X)jeiu s l must be same for all I . Thus, for the 
22 , 55 / : (12,34,13567,24567) configuration, it follows that 

_|_ ^2 _ ^3 _|_ ^4 _ ^1 _|_ ^3 _)_ ^5 _|_ ^6 _|_ ^7 _ ^2 _|_ ^4 _)_ ^5 _|_ ^6 _|_ ^7 


s 3 = s 2 , s 4 = s 1 , s 5 + s 6 + s 7 = 0 , S c = 2(s 1 + s 2 ) 


It can then be shown that s l = s* 

+ m 

satisfy the relations 


C/3* 

CO 

C/3* 

to 

C/3* 

C/3* 

, 5 5 

_l_ _|_ gT _ ~i _|_ ~ 2 

(60) 

Consider equation (34) for s l : 
now gives 

a s l = 

= J2i c d Ri . Applying equations (59) 

a s 1 = Y c 11 Rj 

7 

= c 1 

(R -2 — R' 2 1 + R 5 ~~ R 5 1 ) 


as 2 = J2~ c2Ir i 

7 

= c 1 

(R'2 — Ry — Ro + R 5 ') 


as 3 = 

7 

= c 1 

(—i? 2 + R 2 1 + Ro — R 5 ') 


CK s 4 = E C 47 iU 

7 

= c 1 

(—i? 2 + Ry — Ro + R 5 ') 


as 5 ’ 6 ’ 7 = 1 R i 

= c 1 

(—i? 2 — i? 2 / + R§ + i?5') . 



i 
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The three relations (60) on the s l then imply that (see footnote 11) 

R 

T 


R 2 — R-2' — R5 — R§' — — , R~Y Ri , 

1 


==► s l = 0 ==► s i = S c = J3o = B 1 = 0 . 

Using equations (34) - (36), we then get the same zeroth order results as for 
the (m + 2) dimensional stars where m — 2 now. Namely, we get 


5° = 


R = 


2 lOn (1 - 77) 
1 + UV 


,0 


a = m — 1 + s u = 


m c 

1 + Wv 


as 0 2 m uq (1 — rj) 


(m — 1) e 


2A 0 _ 


= a — c R = 


(1 + W tt ) 2 

V 

(1 + Wir) 2 


where V = (m — 1)(1 + w n ) 2 + 4uq-(l — rj) and m = 2 


Consider the first order equations of motion (31) and (39) - (43). Their 
solutions will give the asymptotic perturbations around the singular solu¬ 
tions. Equations (31) and (57) give 

w n y 1 = — (1 + w n ) u° - cU c + me Y u l ■ 

ie\\i 


Using c 1 = c , c 01 = c° , 
be written as 

ul = 

Zlg. U& 

Zip 


and Bq = £>1 = 0 , the equations for u° and u can 

— Y Ri y 1 + (2S° + m — l) u 

Ul j ; 

8 Y,Ri y 1 + 2 era 

/ 

— Y^y 1 ~ ( m ~ l )u 

m j 


+ aul 


v' 


(c — sc 


ZHi y 1 


(1 + w n ) R 


u° = 


w „ 


— 2 (1 — 77 ) a 


ir 
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The above equation for J2i Ri V 1 follows from Ri — j, the expression for y 1 , 
equation (58), and from the results obtained for R at the zeroth order. 

It is easy to see now that these equations for u° and u are same as those 
for u° and u in the case of the (m + 2) dimensional stars with m — 2 . Hence, 
further analysis of these equations and the consequent results are also the 
same. In particular, u° and u obey equation (46). And, their solutions are 
non oscillatory if the values of the anisotropy parameter q and, equivalently, 
of lie in the ranges given in equation (48). Thus, since m — 2 , certain 
amount of anisotropy, rj > | , is needed to obtain the non oscillatory be¬ 
haviour of the perturbations around the singular solutions in the asymptotic 
region. 

Consider now the equations for u 1 . They are not needed for present 
purposes but we analyse them for the sake of completeness. Since s* = 0 and 
Rj = j, they are given by 

Ka + au's = • ( 61 ) 

4 I 

Evaluating the sum J2i c l1 y 1 using equations (59) gives 12 


E c ^ = 0 =► Y.Pv 1 



For / : (2,2', 5, 5'), the sums J2je\\i u '’ are given by ( u 1 + u 2 ), ( u 3 + w 4 ), 

(u 1 + u 3 + u 5 + u 6 + u 7 ), and {u 2 + u A + u 5 + u 6 + u 7 ) . Using equations (59) 
again gives 


E5 17 V* 

E g3/ y 1 
J2~c AI v 1 
Ec 5 ’ 6 ’ 7 V 


2 c 11 (u 1 — u A ) 

2 c 11 (u 2 — « 3 ) 

2 c 11 (-u 3 — u 2 ) 

2 c 11 (m 4 — u 1 ) 

2 c 11 (u 5 + u 6 + u 7 ) 


(62) 


12 In [37], the factor appearing below was omitted inadvertently. Consequently, the 
equation for (*U given there is incorrect upto this factor in the last term. Equation (63) 
for UU given below is the correct one. 
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It then follows, after some manipulations involving u l s, u l s, equations (61) 
and (62), that 

2 Be 2 

(hW + a(h)a + —(h) = o (63) 

w n 

where (* x ) = ( u 1 — u A ) , (h 2 — u 3 ) , and (u 5 + w 6 + « 7 ) ; and 

(* 2 )^ + c*(* 2)a = 0 (64) 

where (* 2 ) = (u 1 + u 4 ) , (u 2 + u 3 ) , (u 5 + w 6 — 2-u 7 ) , and (h 5 — 2w 6 + u 7 ) . 


7. Conclusion 

We now conclude with a summary and a discussion of some of the issues 
that require further study. 

A brief summary of the present paper is as follows. Unitarity of evolution 
in gravitational collapses implies that horizonless objects must exist which 
can be macroscopic and must be stable. In this paper, with such objects 
in mind, we studied the effects of anisotropy of pressures on the stability of 
stars. The stars are assumed to be static and spherically symmetric in the 
non compact space, to have suitable isometries along the compact directions, 
and to be made up of constituents with linear equations of state. Studying the 
singular solutions and asymptotic perturbations around them, we obtained 
the criteria for the perturbations to be non oscillatory. 

We studied stars in four or higher dimensional spacetime with no compact 
directions, and also two examples of stars in M theory made up of stacks of 
(intersecting) two branes and five branes. A variety of other examples may 
also be studied using the present formulation. We find that non oscillatory 
perturbations around the singular solutions are possible if an appropriate 
amount of anisotropy is present. The details are given in the paper. 

The behaviour of these perturbations lead to corresponding asymptotic 
behaviour in the mass - central density profile of a star enclosed in a box of 
radius r* . The non oscillatory behaviour of the perturbations are likely to 
indicate stability. In that case, singular solutions will correspond to stable 
configurations, and give the mass - radius relation Af(r*) ~ , in the 

limit of large central density or large radius. 

Our results suggest that it may be possible to have stable horizonless 
objects in four or any higher dimensions, and that anisotropic pressures may 
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play a crucial role in ensuring their stability. Although much remains to be 
done, it is worth emphasising that these are important results because they 
bear on the horizonless objects which are implied by unitarity in lien of black 
holes, and they point ont a necessary ingredient for their stability. To actually 
construct such objects, however, requires detailed understanding of many 
issues such as the nature of the constituents and the physical mechanisms 
that may provide the required amount of anisotropy. We now discuss some 
of these issues which may be studied further. 

To show the stability of the equilibrium configurations given by the mass 
- central density profile whose asymptotic behaviour is monotonic non oscil¬ 
latory, one may follow Chavanis and show that the equilibrium configurations 
in the mass - central density profile become unstable beyond the first maxi¬ 
mum; then construct the entire profile and show that it remains monotonic 
and increasing in both the non asymptotic and asymptotic regions. Although 
desireable, we are unable do any of this since the detailed properties of the 
constituents are not known which cause the required amount of anisotropy. 
Proving these things may give valueable insights into the horizonless objects. 

One may try to use scalar fields as in [38] to produce anisotropy and 
thereby to construct a stable horizonless object. The stars constructed in 
these works are unstable for a sufficiently high mass. This may be because 
the scalar field potentials are not tailored to generate anisotropic linear equa¬ 
tions of state. In the cosmological context, a linear equation of state can be 
mimicked using a scalar field with an exponential potential. One may simi¬ 
larly try to mimic anisotropic linear equations of state with scalar fields with 
appropriate potentials and then study the resulting stars. 

There is a vast body of works devoted to construction of anisotropic stars. 
A small sample of them is given in [66] - [73]. Typically, in these works, it is 
found that anisotropy affects stability properties, and that instability sets in 
for a sufficiently high mass. The nature of their ansatzes for anisotropy is very 
different from ours and, hence, there seems to be no discernible contradiction 
between their results and ours. 

There is, however, a distinct possibility that anisotropy of pressures as 
found here is a necessary condition for stability but it may not be suffi¬ 
cient; other ingredient (s) may also be needed. There are two reasons for 
entertaining this possibility. First, in the works [19] - [24] on horizonlcss ob¬ 
jects, anisotropy of the pressures was found to be an important ingredient; 
but a positive cosmological constant, more generally matter with negative 
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pressures, residing in the inner region was also required in an essential way. 
Second, if there is a similarity between the singularities in gravitational col¬ 
lapses and in cosmological big bang/crunch evolutions then the mechanisms 
resolving these singularities may also be expected to be similar. Usually, such 
mechanisms involve new types of matter and/or interactions. For example, 
matter violating null energy conditions can cause a bounce and resolve big 
bang/crunch singularity. Then, going by the similarities, one may also ex¬ 
pect similar ingredients to play a role in stabilising gravitational collapses. 
Thus it seems possible that other ingredient (s), besides anisotropy, may still 
be needed for constructing stable horizonless objects. Nevertheless, using the 
anisotropy criteria given here, one may try to construct such objects and see 
if, and which, further ingredients are needed. 

At a technical level, although our formulations included multi component 
fluids, we only considered cases where Af = 1 , or chose zuf = w n and w 1 = w 
for all / when J\f = 4 . It may be worthwhile to investigate situations 
where more than one component become crucial and play a significant role. 
For example, is it possible that one component dominates the inner regions 
and another the outer regions but such that, together, they lead to stability 
against collapse? 
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